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Abstract. - We study the orientational dynamics of infinitely thin hard rods of length L, with 
the centers-of-mass fixed on a simple cubic lattice with lattice constant a. We approximate the 
influence of the surrounding rods onto dynamics of a pair of rods by introducing an effective 
rotational diffusion constant D(l), I = L/a. We get D(l) oc [1 — v(l)], where v(l) is given 
through an integral of a time-dependent torque-torque correlator of an isolated pair of rods. 
A glass transition occurs at l c , if v(l c ) = 1. We present a variational and a numerically exact 
evaluation of v(l). Close to l c the diffusion constant decreases as D(l) oc (Z c — l) 1 , with 7 = 1. 
Our approach predicts a glass transition in the absence of any static correlations, in contrast 
to present form of mode coupling theory. 



When a system of interacting species is cooled sufficiently fast it will typically undergo a 
glass transition. In the last two decades there has been a significant progress in microscopic 
understanding of this transition. The first approach is the so-called mode coupling theory 
(MCT) P| , first suggested and mainly worked out by Gotze and coworkers. For reviews 
see Refs. |3|, [I|,f|. MCT yields a closed set of equations for time-dependent correlators. This 
set involves as an input static correlators which depend only on thermodynamic variables like 
temperature T, density n, etc. With decreasing T (or increasing n) local ordering described 
by peaks of the static correlators grows. When the peaks' magnitude becomes large enough 
MCT predicts an ergodicity breaking transition at a critical temperature T c (or density n c ). 
This dynamic transition is interpreted as an ideal glass transition. 

A different approach has recently been proposed by Mezard and Parisi HQ]. Inspired by 
the replica theory for spin glasses they developed a microscopic replica theory for structural 
glasses. This theory predicts a static glass transition at a temperature Tf which is lower than 
T c . One of its characteristic features of this transition is the vanishing of the configurational 
entropy at Tf. The physical picture behind that replica theory was suggested earlier by Singh 
et al. f§. 

Despite providing new valuable insights, these microscopic theories can not be complete, 
because there exists a purely dynamical mechanism which leads to a glass transition which is 
not based on the existence of any static correlations. As an example, let us consider a simple 
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cubic lattice with N lattice sites and lattice constant a with hard rods of length L and diameter 
d with their centers fixed on the lattice sites (Fig. [l]). This model system was first introduced 
by Rubinstein and Obukhov B. It can be used to describe orientational glasses pi 111]] F) 



The increase of the reduced rod length I = L/a leads to an increase of the steric hindrance, 
which may lead to a glass transition. Now suppose that the thickness d goes to zero, for fixed 
a and L. For d = all static orientational correlations vanish and there is no static phase 
transition. Nevertheless, the simulational studies by Renner et al. |l3| for a fee lattice and 
Obukhov et al. j[4| clearly demonstrate the existence of a (dynamical) glass transition at a 
critical length l c = L c /a. In addition, it is particularly interesting that the time-dependent 
orientational correlators studied in Ref. [ p"3| exhibit a two-step relaxation process similar to 
that observed for supercooled liquids and predicted by MCT [||,|],|5|]. However, as Renner et 
al. [fl3| stressed the MCT in its present form fails to describe this dynamical transition. The 
reason is obvious: the vertices which enter into the MCT equations vanish pj . Thus, present 
formulation of MCT as well as the replica theory cannot explain the numerical findings 
of Ref. Jl^,|lJ], due to the complete trivial statics for d = 0. 

This discussion demonstrates the necessity to derive a microscopic theory which is capable 
to describe a purely dynamical glass transition for systems with completely trivial statics. 
To achieve this for the infinitely thin hard rods on a cubic lattice is the goal of the present 
contribution. Details and additional features of our approach will be given elsewhere ]l6[ . 

Let us consider the time-dependent 1-rod density pn} (tyt) for a tagged rod at lattice site 
n\ with orientation f2 = (O, <f) and 

P W(n;o) = W) (i) 

as initial condition, where 5(f2|0') = sin 0(5(9 — &)S(4> — (/>'). If / < 1 the rods cannot 
touch and the tagged rod will perform a free rotational motion. However, for I > 1 steric 
hindrance exists and the tagged rod's motion becomes influenced by its nearest, next-nearest, 
etc. neighbours. These interactions will influence the rotational motion of the rod at site n\. 
With increasing I rotational diffusion constants will decrease and finally may become zero at 
a critical length l c . 

The time derivative of pni(Q',t) is related to the 1-rod current density j!V(^;i) by the 
continuity equation: 

|pW(fi;0 + Vo-£ 1 ) («;i) = 0. (2) 
We define a generalized rotational diffusion tensor, D(f2, Q';t), by: 



j (1 



l ](n ; t) = -J dt' J dn' D(n,fi' ; t-f) • v^/iVC^';*')- (3) 



In the following we derive a closed equation for the diffusion tensor D which will allow 
to investigate its ^-dependence. To derive this equation we use Smoluchowski (Brownian) 
dynamics (glass transition should not depend on the microscopic dynamics [|~7], E0§)- We 
start from the so-called generalized Smoluchowski equation for the iV-rod probability density 

( 1 ) We stress that this type of model is quite different to the mixed crystals studies by Michel |jl2t where the 
glass transition occurs due to quenched disorder 

( 2 ) There exist non-generic model systems for which the vertices do not vanish although there are no static 
correlations |15|. 
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p N (n lr --,n N ;ty. 
d 



N 



dt 



p N (n 1: ---,n N ;t) = D J2 {ni,---,^N)]PN(flu--'^N',t) (4) 



n=l 



where Dq is the bare rotational diffusion constant, V ni = Vn n , and 

T n (fti, • ■ • , ftjv) = ^ Tj n (flj,rt n ) 



where 



T y -(ft,, ft,-) = % ( Ui x f+) 5(r+ - 0+)6 - - \s 



L 



(5) 



(6) 



describes the singular torque when two rods at sites i and j with vector distance touch 
each other [jl9| Q . is the component of perpendicular to the plane defined by the unit 
vectors along the rods, Ui and Uj, = |r^ | and r i:) = r^/r^ . Sjj and Sji (which depend on 
fti, ft,) is the distance of the contact point to the center of rod i and rod j, respectively (cf. 
Fig. |^). S(x — + ) stands for lim S(x — d). The reader should note that is t-independent 

and is equal to the lattice vector = R; — Rj, in contrast to a fluid of hard rods (needles). 

A generalized Smoluchowski equation similar to Eq. (3) had been used to describe a fluid 
of infinitely thin rods with randomly frozen orientations [2y,^T|. In this case there is no glass 
transition if the longitudinal component of the bare translational diffusion tensor is positive. 

From Eq. (||) one can derive a hierarchy of equations of motion for the reduced j-rod 
densities /»^l.. n . (fti, ■ ■ ■ , ft.,-; t) p2fl. For instance, for j = 1 one gets: 



-p«(ft 1 ;t)= J DoV ni 



V„ l/9 W(ft i; ;)-E / dft 2 T„ iri2 (ft 1 ,ft 2 )p( l 2 i )„ 2 (ft 1 ,ft 2 ;0 



(7) 



The first term on the r.h.s. describes the free Brownian motion. The second term describes 
the influence of a rod at site n 2 on the tagged rod. It involves the 2-rod density p ni n 2 (ftii ft 2 ; t). 
To proceed it is convenient to introduce the fluctuations: 



6p 



&?... nj («!,..., n i5 t) = P { il. n] (ftx, • • • , ft i; t) - jAjzrg^l-n, (n Xj ■ • • , n 3 ) P ^(n i; t) (8) 



where gn}... nj (fti, • ■ ■ , ftj) is the equilibrium j-rod distribution function. Note that the equi- 
librium distribution functions are equal to one almost everywhere. 

Using Eq. (ph one obtains a hierarchy of equations for 5pV>. The lowest level of this 
hierarchy is still given by Eq. (Q) where one is allowed to replace p£}n 2 by 8pn}n 2 , because 
Yj J di} 2 T ni „ 2 (fti,ft 2 ) gnin 2 (tti, ft 2 ) = 0. At the second level one gets @: 



^ ^„ 2 (ft! , ft 2 ; t) = - i- (V ni ^ (ftx , ft 2 )) • j « (fti ; t) 



i>o {v ni • [v ni - T„ 1Il2 (fti,ft 2 )] + (i «-> 2)} <^„ 2 ^ ,n 2; t) + ^w na (n 1 ,n 2 ;t) (9) 



( 3 )Since V u „ is used in Ref. Il9| instead of Vn n their expression for Tij differs slightly from ours. 
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where 



-d J2 v ^-12 f dn 3 T n ^ 3 (n v ,n 3 )s P ^ n2n3 (n 1: n2,n 3 ;t). (10) 

Ailing describes the influence of a third rod at n 3 on the rods at n\ and n 2 . Here we follow 
Rcfs. [^o[|l[ and approximate this influence by introducing an effective diffusion tensor D eff . 
Specifically, in Eq. (|J) we replace Dq by the effective diffusion tensor and we neglect An}m '■ 

{A>v ni • [v ni - T nina (ni,n a )] + (i <-> 2)}6p^ n2 (n 1 ,n 2 ;t) -> (ft) 

v ni - J dt' j dn[D eS (n 1 ,n[-t-t')-[w nti -T nin2 (n[,n 2 )] + ^ 
^(ni.najtjwo. (12) 

In order to get a self consistent equation we then set in Eq. (O) D cff = D. This approximation 
leads to a closed set of equations, given by Eqs. (||), (ft), (|ll|) and (^2|). Taking Laplace 
transforms of these equations one finally gets a self consistent equation for the diffusion tensor 
D(ili,£7 2 ) = lim D(f2i,f2 2 ; z), where the latter is the Laplace transform of D(f2i,f2 2 ;£) flq| . 

z— >0 — 

For technical convenience we perform here an additional approximation: 

D^coi.n;) «D(i)<f a/s *(n 1 |ni). (13) 

Substitution of Eq. (|l3|) into the self consistent equation for D yields Jl6| : 

D(l)=D [l-v(l)] (14) 
with the /-dependent coupling function 

oo 

v(l) = lY\ / *(T nin2 (0 1 ,Q 2 )-e £(2)tt T nin2 (n 1 ,fi 2 )) . (15) 
where £( 2 ^ is an adjoint Smoluchowski operator for an isolated pair of rods 

£ (2)t = [v ni +T ni „ 2 (ni,n 2 )] • v ni + [v„ a +T„ 2ni (n 2 ,ni)] • v„ 2 . (16) 

The time-dependent torque-torque correlator on the r.h.s. of Eq. ( |l5| ) is defined by: 
(T nin2 (Ox,0 2 ) • e £<2)tt T nin2 (nx,0 2 )) 

= jkf J dfl1 /^P^L^i.^T^^.O^-e^'^T^^!,^). (17) 

Due to the lattice translation invariance v(l) does not depend on n\. Moreover, v(l) > 0, 
where the equality sign holds only for I < 1 and v(l) ~ I 3 for 2 ;» 1. 



R. SCHILLING and G. SZAMEL: MICROSCOPIC THEORY FOR THE GLASS TRANSITION IN A SYSTEM WITHOUT STATIC CORRELATIONS^ 



The coupling constant v(l) can be calculated via a Brownian dynamics simulation of 
the isolated two-rod system or it can be estimated analytically. The former approach is 
straightforward. In the following we briefly discuss the latter. Eq. ( |T5| ) can be rewritten as 

v(i) = ^(T„ in2 (o 1 ,o 2 ) • (-rWtj-iTn^cni.na)). (is) 

"2 

We introduce a vector function f ni „ 2 (fii, fi 2 ) such that 

C m f nin2 (0 1 ,0 2 ) = T nin2 (f2 1 ,0 2 ) . (19) 
If f is a solution of Eq. ([l9|), then v(l) can be readily calculated: v(l) = — ^ 5Z(T ni „ 2 (£~2i , O2) ■ 

"2 

f ni „ 2 (f2i, f2 2 )}- fnin 2 (^lj ^2) can be determined from the variational equation of the following 
functional (for a similar discussion see Ref. p3|): 

^[{^(Q^O,)}] 

= ^{f nin2 (Cl u n 2 ) ■ C^fn^Q!,^)) - ?(f ni „ 2 (Q l5 n 2 ) ■T niva (Qx,^h)}- (20) 

v(l) is the maximum value of T and v™(l) = ^"[{f™^ 2 (^1, ^2)}] < f (0 for all trial functions 
f«in 2 ■ We have used the following simple trial function 

(cos9icos0i \ f , , . , , 

Sill KJ ^ / 

where A is a variational parameter. Note that f™^ 2 is independent on m, n 2 and does not 
depend explicitly on I. For this simple trial function the ^-dependence of v vav (l) comes from 
that of T„ l7l2 

In Fig. ||we present v(l) for a simple cubic lattice as obtained from Eq. ( |l5| ) by a computer 
simulation of the isolated 2-rod system and from the variational approach. The glass transition 
is defined by D(l c ) = which results in v(l c ) — 1. A finite value for l c > 1 must exist because 
v(l) > and v(l) ~ I 3 for I > 1. From Fig. | we get l™ m 3.45 and I™ ^ 5.7, respectively. 
For the fee lattice we have found |§ Z™ 1 ^ 2.2 and l v c aT S 3.96 which reasonably fits S 2.7 
|( 4 )| from the MD-simulation p3[ . From this point of view our microscopic approach, including 
all approximations, leads to a satisfactory value for the critical length. Ref. jl3| has also 
demonstrated that the rotational diffusion constant follows a power law D(l) ~ (l c — Z) 7 , 
close to Z c with exponent 7 MD = 4.2. Our approach yields a power law with 7 = 1, which 
significantly deviates from the MD result. This discrepancy may be due to some uncertainties 
of the simulation (7 MD = 4.2 is an unusually high value; for supercooled liquids 7 is close to 
two) or it may originate from our approximations. With regard to the latter it should be 
noted that the simplest version of Edwards and Vilgis' description of the glass transition in 
a dense fluid of rods of finite thickness |24j predicts a power law decay of the self-diffusion 
coefficient with 7=1. This result is changed in more complicated formulations of the theory. 
Whether our approach can be developed in an analogous way is left for the future study. 

To summarize, we have considered a system of hard rods fixed with their centers on a simple 
cubic lattice. We approximated the influence of a third rod on a pair of rods by introducing 



( 4 )Obukhov et al. considered a system of rods attached to the lattice sites of an sc lattice at one of their 
end points and obtained a slightly different value, l c = 4.5. 
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an effective diffusion tensor and derived a closed set of equations for the 1- and 2-rod density. 
These equations lead to a self consistent equation for a rotational diffusion constant D(l). In 
contrast to MCT it is not a static but a dynamical 2-particle correlator which determines the 
slowing down of the dynamics and finally the glass transition. Nevertheless, similar to the 
MCT, a cage effect exists in our model, too, but the cage is of pure dynamical nature. If I is 
large enough, the tagged rod will perform librations within this "dynamical cage" built by all 
the neighbouring rods, without existence of any static correlations. Note that this mechanism 
exists even for rods with finite d, as long as the static correlations are weak enough. With 
increasing d there will be a crossover to the "static" cage effect as described by MCT. Whether 
or not MCT can be extended such that this "dynamical" cage effect is included is, of course, 
a challenge for the future. In any case, the present model exhibits glassy behavior without 
possessing any static correlations. 
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Fig. 1 - Unit cell of the cubic lattice with hard rods of length L and diameter d with their centers at 
the lattice sites. 




Fig. 2 - Illustration of the geometrical quantities defined in the text for two needles i and j. 
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Fig. 3 - v(l) for a simple cubic lattice from the numerically exact (squares) and the variational 
calculation (solid line). l c denotes the critical length for which v(l c ) = 1. 



